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Connecting Polyakov Loops to Hadrons 
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The order parameter for the pure Yang-Mills phase transition is the Polyakov loop, which encodes the symmetries of 
the Zn center of the SU (N) gauge group. The physical degrees of freedom of any asymptotically free gauge theory are 
hadronic states. Using the Yang-Mills trace anomaly and the exact Zn symmetry we show that the transfer of information 
from the order parameter to hadrons is complete. 



1 Introduction 

Recently in we have analyzed the problem 

of how, and to what extent the information encoded in 
the order parameter of a generic theory is transferred 
to the non-order parameter fields. This is a fundamen- 
tal problem since in nature most fields are non-order 
parameter fields. 

This problem is especially relevant in QCD and QCD- 
like theories since there is no physical observable for 
deconfinement which is directly linked to the order 
parameter field. Here we mainly review the first 
work [in on this issue, which deals with the confine- 
ment/deconfinement phase transition in pure Yang- 
Mills theory. 

Importance sampling lattice simulations are able to 
provide vital information about the nature of the tem- 
perature driven phase transition for 2 and 3 color 
Yang-Mills theories with and without matter fields 
(see [ El for 3 colors) . At zero temperature the 
SU{N) Yang-Mills theory is asymptotically free, and 
the physical spectrum of the theory consists of a tower 
of hadronic states referred to as glueballs and pseudo- 
scalar glueballs. The theory develops a mass gap and 
the lightest glueball has a mass of the order of few 
times the confining scale. The classical theory pos- 
sesses conformal symmetry, while quantum corrections 
lead to a non- vanishing trace of the energy momentum 
tensor. 

At nonzero temperature the Zn center of SU (N) is a 
relevant global symmetry [[7], and it is possible to con- 
struct a number of gauge invariant operators charged 
under Z^- Among these the most notable one is the 
Polyakov loop: 
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V denotes path ordering, g is the SU {N) coupling con- 
stant, X is the coordinate for the three spatial dimen- 
sions, and r is the Euclidean time. The £ field is real 
for N = 2, while otherwise complex. This object is 
charged with respect to the center Zjy of the SU{N) 
gauge group [|7], under which it transforms as £ ^ zl 
with z G Zn. A relevant feature of the Polyakov loop 
is that its expectation value vanishes in the low tem- 
perature regime, and is non-zero in the high temper- 
ature phase. The Polyakov loop is thus a suitable or- 
der parameter for the Yang-Mills temperature driven 
phase transition [[7]. 

Here we consider pure gluon dynamics. This allows us 
to have a well defined framework where the Zn sym- 
metry is exact. The hadronic states are the glueball 
fields (H) and their effective theory at the tree level is 
constrained by the Yang-Mills trace anomaly. 

A puzzle is how the information about the Yang- 
Mills phase transition encoded, for example, in the 
Zn global symmetry can be communicated to the 
hadronic states of the theory. Here we propose a con- 
crete way to resolve the puzzle. 

As basic ingredients we use the trace anomaly and the 
Zn symmetry. They will be enough to demonstrate 
that the information carried by £ is efficiently and com- 
pletely transferred to the glueballs. More generally, 
the glueball field is a function of £: 

H^H[£]. (1) 

Our results can be tested via first principle lattice sim- 
ulations. 

2 The Basic Properties 

The hadronic states of the Yang-Mills theory are the 
glueballs. At zero temperature the Yang-Mills trace 
anomaly has been used to constrain the potential of 
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the lightest gluebaU state [El: 



mi = fin 



H 

A4 



(2) 



A is chosen to be the confining scale of the theory 
and H is a mass dimension four field. This poten- 
tial correctly saturates the trace anomaly when H is 
assumed to be proportional to Tr[GpyG^''] and G^i/ 
is the standard Yang-Mills field strength. The poten- 
tial nicely encodes the properties of the Yang-Mills 
vacuum at zero temperature and it has been used to 
deduce a number of phenomenological results Ef- 
fective Lagrangians play a relevant role for describing 
strong dynamics. Recently, for example a number of 
fundamental results at zero temperature and matter 
density about the vacuum properties and spectrum of 
QCD have been uncovered [[TU]. 

At high temperature the Yang-Mills pressure can 
be written in terms of the field i according to the 
Polyakov Loop Model (PLM) of [Ej. This free energy 
must be invariant under Zn and it takes the general 
form: 



V[£] ^T'^T[l]. 



(3) 



is a polynomial in ^ invariant under Z^r, and its 
coefficients depend on the temperature, allowing for a 
mean field description of the Yang- Mills phase transi- 
tions. 

We now marry the two models by requiring both fields 
to be present simultaneously at non zero tempera- 
ture. The theory must reproduce the ordinary glue- 
ball Lagrangian at low temperatures and the PLM La- 
grangian at high temperatures. In [ 1 the following 
effective potential was proposed: 



y[i/,£] = f In 



H 

A4 



Vt[H] + HV[l\+T^V [^],(4) 



where V [t] and V [I] are general (but real) polynomi- 
als in £ invariant under whose coefficients depend 
on the temperature. The explicit dependence is not 
known and should be fit to lattice data. Dimensional 
analysis and analyticity in H when coupling it with 
£ severely restricts the effective potential terms. We 
stress that HV[l] is the most general interaction term 
which can be constructed without spoiling the zero 
temperature trace anomaly. Further nonanalytic in- 
teraction terms can arise when considering thermal 
and quantum corrections. These are partially con- 
tained in Vt[H] which schematically represents the 
temperature of a gas of glueballs. In the following 
we will not investigate in detail such a term. Our the- 
ory cannot be the full story since we neglected (as cus- 
tomary) all of the tower of glueballs and pseudo-scalar 



glueballs as well as the infinite series of dimensionless 
gauge invariant operators with different charges with 
respect to Z^. Nevertheless, the potential is suffi- 
ciently general to capture the essential features of the 
Yang-Mills phase transition. 

When the temperature T is much less than the confin- 
ing scale A the last term in Eq. Q can be safely ne- 
glected. Since the glueballs are relatively heavy com- 
pared to the A scale their temperature contribution 
Vrt-ff] can also be disregarded. At low temperatures 
the theory reduces to the standard glueball potential 
augmented by the third term which does not affect the 
trace anomaly. 

At very high temperatures (compared to A) the last 
term dominates {H itself is very small) recovering the 
picture in which t dominates the free energy. In this 
regime we have J-[P\ = V[f]- 

We can, in principle, compute all the relevant ther- 
modynamic quantities in our approach, i.e. entropy, 
pressure etc., but this is not the main scope of this 
work. 

A relevant object is the trace of the energy-momentum 
tensor At zero temperature, and when the poten- 
tial is a general function of a set of bosonic fields 
with mass-dimensions (i„, one can construct the asso- 
ciated trace of the energy-momentum tensor via: 



% = 4F[$„] 



E 



(5) 



At finite temperature we still define our temperature 
dependent energy-momentum tensor as in Eq. |(SJ). 
Here H possesses engineering mass dimensions 4 while 
i is dimensionless, yielding the following temperature 
dependent stress energy tensor: 



e^(T) 



-2H + AT'^V[ 



1 - H- 



SH 



Vt[H]. (6) 



is normalized such that (O|0^|O) — e ~ 3p with 
e the vacuum energy density and p the pressure. At 
zero temperature only the first term survives, yield- 
ing magnetic type condensation typical of a confining 
phase, while at extremely high temperature the sec- 
ond term dominates, displaying an energy density and 
a pressure typical of a deconfined phase. 

The theory containing just £ can be obtained integrat- 
ing out H via the equation of motion: 



SV[H, £] 
SH 



0. 



(7) 



Formally this is justifiable if the glueball degrees of 
freedom are very heavy. For simplicity we neglect the 
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contribution of Vt[^^], as well as the mean- field the- 
ory corrections for i. However, a more careful treat- 
ment which also includes the kinetic terms should be 
considered Within these approximations the 

equation of motion yields: 



A* 

H\l] = — exp[-2P[^]] 
e 



(8) 



The previous expression shows the intimate relation 
between t and the physical states of strongly interact- 
ing theories. 

After substituting Eq. (jHJ back into the potential Q 
and having neglected Vri-ff] we have: 



V[l] = T^V[ 



A4 
2e 



exp[-2P[^]] 



(9) 



This formula shows that for large temperatures the 
only relevant energy scale is T and we recover the PLM 
model. At low temperatures the scale A allows for new 
terms in the Lagrangian. Besides the and the A^ 
terms we also expect terms with coefficients of the type 
TK^ and T^A^ and T'^A. However in our simple model 
these terms do not seem to emerge. 



Expanding the exponential we have: 



V[i] = T^V[l] + —V[l] 



A4 
2e 



(10) 



Since V[P\ and 'P[P\ are real polynomials in d invariant 
under we immediately recover a general potential 
in I. 



3 The two Color Theory 

To illustrate how our formalism works we first consider 
in detail the case N = 2 and neglect for simplicity the 
term Vr [H] . This theory has been extensively studied 
via lattice simulations [ I12L [T5| and it constitutes the 
natural playground to test our model. Here € is a real 
field and the Z2 invariant V[£] and V[£] are taken to 
be: 



v[e] 
r\£] 



air 



- a2£^ + 0(^^ 
Oil'), 



(11) 



with fli, 02 and 61 unknown temperature dependent 
functions, which should be derived directly from the 
underlying theory. Lattice simulations can, in princi- 
ple, fix all of the coefficients. In order for us to investi- 
gate in some more detail the features of our potential, 
and inspired by the PLM model mean-field type of 



approximation, we first assume 02 and bi to be posi- 
tive and temperature independent constants, while we 
model oi = a(r* — T)/T, with a constant and 
a another positive constant. We will soon sec that 
due to the interplay between the hadronic states and 
£, need not to be the critical Yang-Mills tempera- 
ture while ai displays the typical behavior of the mass 
square term related to a second order type of phase 
transition. 

The extrema are obtained by differentiating the po- 
tential with respect to H and £: 



dV 

dV_ 
'dl 



In 
2" 



eH 



A4 



H 



eH 



A4 



2£T^ [ ai -t- ^61 + 2a2£^ 



f = 0(12) 
0. (13) 



3.1 Small and Intermediate Temperatures 

At small temperatures the second term in Eq. (|13|) 
dominates and the only solution is £ — Q. A vanishing 
£ leads to a null 'P[£] yielding the expected minimum 
for H: 



A* 
e 



(14) 



Here £ and H decouple. 



We now study the solution near the critical tempera- 
ture for the deconfinement transition. At all temper- 
atures for which 



T^ai + Hbi = T^a{T^ - T) + Hbi > 0, 



(15) 



the solution for £ is still £ = 0, yielding Eq. lO. The 
critical temperature is reached for 



b^A^ 

ea ' 



(16) 



and can be determined via lattice simulations. We see 
that within our framework the latter is related to the 
glueball (gluon-condensate) coupling to two Polyakov 
loops and it would be relevant to measure it on the 
lattice. At T = Tc, ^ = and H = K'^/e. 



Let us now consider the case T = Tc + AT with 
AT 



Tr 



< 1. 



(17) 



Expanding (^)^ at the leading order in l\T jT^ yields: 
l + 3|i|i AT 



{£Y^ 



'2a2 1 _ ilAl Tc 

ea2 Ti 



(18) 
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We used Eq. lfTH|l and Eq. ifT^ which relates the tem- 
perature dependence of H to the one of i. At high 
temperatures (see next subsection) {£) can be normal- 
ized to one by imposing ajlai — 1 and the previous 
expression reads: 



1 



AT 

ea 



4Tc - 



(1 - 26i)re 



AT 



(19) 



For a given critical temperature consistency requires 
h\ and T* to be such that: 



4Tc - ST* 



(1 - 26i)rc + 26ir* 



> 0. 



(20) 



In this regime, the temperature dependence of the 
gluon condensate is: 



A-* 

i?) = — cxp [-26i(^)2] 



(21) 



We find the mean field exponent for i.e. f ^ increases 
linearly with the temperature near the phase transi- 
tion ^. Interestingly, the gluon-condensates drops ex- 
ponentially. This drop is triggered by the rise of I and 
it happens in our simple model exactly at the deconfin- 
ing critical temperature. Although the drop might be 
sharp it is continuous in temperature. This is related 
to the fact that the phase transition is second order. 
Our findings strongly support the common picture ac- 
cording to which the drop of the gluon condensate 
signals, in absence of quarks, deconfinement. 

3.2 High Temperature 

At very high temperatures the second term in Eq. H13|) 
can be neglected and the minimum for I is: 



(22) 




For iJ we have now: 



A" 

(ff ) = — exp 

e 



-26i^ 
2a2 



A4 

= — exp [-2&i] . (23) 

e 



In the last step we normalized {I) to unity at high 
temperature. In order for the previous solutions to be 
valid we need to operate in the following temperature 
regime: 



V a 



(24) 



^Corrections to the mean field are large and must be taken into 
account. 



We find that at sufficiently high temperature (iJ) is 
exponentially suppressed and the suppression rate is 
determined solely by the glueball ~ mixing term 
encoded in V\(\. The coefficient h\ should be large 
(or increase with the temperature) since we expect 
a vanishing gluon-condensate at asymptotically high 
temperatures. Clearly it is crucial to determine all of 
these coefficients via ffi'st principle lattice simulations. 
The qualitative picture which emerges in our analysis 
is summarized in Fig. ^ 




0.98 0.99 



Figure 1. The thin line is the gluon condensate {//) nor- 
malized to A*/e as function of temperature. The thick line 
is {€) normalized to unity as function of temperature. We 
chose for illustration a = 1, 6i = 1.45 and Tc ~ 1.16A. 



4 The three color theory 

Zz is the global symmetry group for the three color 
case and £ is a complex field. The functions and 
V\l\ are: 



6i|£p - 



-021^1' 



(25) 



with oi, 02, 03 and h\ unknown temperature depen- 
dent coefficients which can be determined using lattice 
data. Here we want to investigate the general relation 
between glueballs and £, so we will not try to find the 
best parameterization to fit the lattice data. In the 
spirit of the mean field theory we take 02, 03 and h\ to 
be positive constants while ai = a(T, — T)/T. With 
I = \P\e'^'^ the extrema are now obtained by differenti- 
ating the potential with respect to \l\ and Lp: 



dV 
dH 
dV 

W\ 

dV 

dip 



In 
"2" 



eH 



A4 



eH 



A4 



+ 6iKP = 



H 



= 2\£\T'' l^ai + —h 



3a3\£\ cos{3ip) + 2a2| 
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At small temperature the H/T'^ term in the second 
equation dominates and the solution is (|-^|) = 0, 
{H) = A'^/e. The last equation is verified for any 
(</?), so we choose {^p) = 0. The second equation can 
have two more solutions: 



3 03 

4 02 



± 



16 a2 



a{T - n) biH 



2Ta. 



2a2r4' 



(27) 



whenever the square root is well defined, i.e. at suf- 
ficiently high temperatures. The negative sign corre- 
sponds to a relative maximum, while the positive one 
to a relative minimum. We have then to evaluate the 
free energy value (i.e. the effective thermal potential) 
at the relative minimum and compare it with the one 
at £ = 0. The temperature value for which the two 
minima have the same free energy is defined as the 
critical temperature and is: 



bi A4 
ear? 



aa2 



aa2 



(28) 



When 03 vanishes we recover the second order type 
critical temperature. To derive the previous expres- 
sion we held the value of H fix to A^/e at the tran- 
sition point. In a more refined treatment one should 
not make such an assumption. Below this tempera- 
ture the minimum is still for {£) = and (H) = A'^/e. 
Just above the critical temperature the fields jump to 
the new values: 



03 
a2 ' 



A 



H) = _exp[-2&i(|£|)2]. 



(29) 



Close to, but above (i.e. T :^Tc + AT) we have: 
i\l\)^^ + p^, (30) 

0-2 J-c 

with 



representing the behavior of the Polyakov loop and the 
gluon condensate is presented in Fig. |21 

Since we are in the presence of a first order phase 
transition higher order terms in Eq. (|25|l may be im- 
portant. Lattice simulations have shown however, that 
the behavior of the Polyakov loop for 3 colors resemble 
a weak first order transition (i.e. small 03) partially 
justifying our simple approach. The approximation 
for our coefficients is too crude and it would certainly 
be relevant to fit them to lattice simulations. 

What we have learned is that the gluon condensate, 
although not a real order parameter, encodes the in- 
formation of the underlying Z3 symmetry. More gen- 
erally, we have shown that once the map between 
hadronic states and the true order parameter is known, 
we can use directly hadronic states to determine when 
the phase transition takes place and also determine 
the order of the phase transition. For example by 
comparing Fig. ^ and Fig. El we immediately notice 
the distinct behaviors in the temperature dependence 
of the gluon condensate near the phase transition. 



0.6 
0.4 
0.2 





<H> 







0.99 



1.01 1.02 1.03 

T/Tc 



Figure 2. A cartoon sketching the gluon condensate {H) 
normalized to A*/e and the (|^[) (thick line) as function of 
temperature. We chose for illustration as = 0.3, 02 = 1, 
a = 2, 61 = 0.7 and Tc = 1.2A. 



aa2 



AnTr — 3T* 



as a'2Tc - Ahia{nTc - T,) ' 
aa2 -\- al 



aa2 



(31) 



a positive function of the coefficients of the effective 
potential. In this regime 



A* 

(H) = — exp 

e 



„, ,03 AT, 9 
a2 ic 



(32) 



At high temperature we expect a behavior similar to 
the one presented for the two color theory. A cartoon 



5 Conclusions 

Our new theory is able to account for many features 
inherent to the Yang-Mills deconfining phase transi- 
tion. We related two very distinct and relevant sec- 
tors of the theory: the hadronic sector (the glueballs), 
and some dimensionless fields {€} charged under the 
discrete group Zjq understood as the center of the un- 
derlying SU{N) Yang-Mills theory. 

The gluon-condensate is, strictly speaking, not an or- 
der parameter for the deconfining Yang-Mills phase 
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transition. However we have shown that the infor- 
mation encoded in the true order parameter t is effi- 
ciently communicated to the gluon condensate. Since 
the exponential drop of the condensate just above 
the Yang-Mills critical temperature is a direct con- 
sequence of the behavior of the true order parameter 
at the transition we can consider this drop as a strong 
signal of deconfinement. This drop has already been 
used in various models for the Yang-Mills phase tran- 
sition. We have also seen that the reduction in the 
gluon-condensate is associated to the increase of the 
Polyakov loop condensate I. The information about 
the order of the phase transition is also transferred 
to the behavior of the gluon condensate. Indeed, from 
Fig. Hand Fig.[21we deduce that the drop is continuous 
for the gluon condensate in the two color case, while is 
discontinuous for the three color theory. We now have 
a deeper understanding of the mechanism for transfer- 
ring information from the Yang-Mills order parameter 
to the physical states. Other theoretical investigations 
also seem to support the present results [I14|. 

It is worth mentioning that the Polyakov loop need 
not to be the only acceptable order parameter. For 
example using an abelian projection one can define a 
new (non local in the cromomagnetic variables) order 
parameter [ I15j . Our model can be, in principle, mod- 
ified to be able to couple the hadronic states to any 
reasonable Yang-Mills order parameter. 

Finally the present approach has also been extended 
to strongly interacting theories with fermions in the 
fundamental or adjoint representation of the gauge 
group. Within our framework we were then able to 
understand the relation between confinement and chi- 
ral symmetry breaking, and to make a number of pre- 
dictions which can be tested via first principle lattice 
simulations [0]. 

It is now clear, that the critical behavior defined and 
governed by the order parameter can be studied via 
the singlet fields. 
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